APS preprint 



Distribution of the Largest Aftershocks in Branching Models of 
Triggered Seismicity: Theory of Bath's law 

A. Saichev^'^ and D. Sornette^'^ 

^Mathematical Department, Nizhny Novgorod State University, 
Gagarin prosp. 23, Nizhny Novgorod, 603950, Russia 
^Institute of Geophysics and Planetary Physics, 
University of California, Los Angeles, GA 90095 
^Institute of Geophysics and Planetary Physics and Department of Earth and Space Sciences, 

University of Galifornia, Los Angeles, GA 90095 
^Laboratoire de Physique de la Matiere Gondensee, 
GNRS UMR 6622 and Universite de Nice-Sophia Antipolis, 06108 Nice Gedex 2, FrancJ*| 

(Dated: February 2, 2008) 

Abstract 

Using the ETAS branching model of triggered seismicity, we apply the formalism of generating prob- 
ability functions to calculate exactly the average difference between the magnitude of a mainshock and 
the magnitude of its largest aftershock over all generations. This average magnitude difference is found 
empirically to be independent of the mainshock magnitude and equal to 1.2, a universal behavior known 
as Bath's law. Our theory shows that Bath's law holds only sufficiently close to the critical regime of 
the ETAS branching process. Allowing for error bars ±0.1 for Bath's constant value around 1.2, our 
exact analytical treatment of Bath's law provides new constraints on the productivity exponent a and 
the branching ratio n: 0.9 ^ a < 1 and 0.8 ^ n < 1. We propose a novel method for measuring a 
based on the predicted renormalization of the Gutenberg-Richter distribution of the magnitudes of the 
largest aftershock. We also introduce the "second Bath's law for foreshocks: the probability that a main 
earthquake turns out to be the foreshock does not depend on its magnitude p. 
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I. INTRODUCTION 



This paper is part of our continuing effort to develop a complete theory of seismicity within 
models of triggered seismicity, which allows one to make quantitative predictions of observables 
that can be compared with empirical data . We study the general branching process, 

called the Epidemie-TVpe Aftershock Sequence (ETAS) model of triggered seismicity, introduced 
by Ogata in the present form [7] and by Kagan and Knopoff in a slightly different form [S| and 
whose main statistical properties are reviewed in j^. The ETAS model belongs to a general class 
of branching processes ^, Q| , and has in addition the property that the variance of the number of 
earthquake progenies triggered in direct lineage from a given mother earthquake is mathematically 
infinite. This model has been shown to constitute a powerful null hypothesis to test against 
other models Q|. The advantage of the ETAS model is its conceptual simplicity based on three 
independent well-found empirical laws (Gutenberg-Richter distribution of earthquake magnitudes, 
Omori law of aftershocks and productivity law) and its power of explanation of other empirical 
observations (see for instance p| and references therein). 

Here, we develop a theoretical formulation based on generating probability functions (GPF) to 
construct the distribution of magnitudes of the largest triggered event (largest aftershock) within 
the cascade comprising all triggered events of a given source earthquake. This allows us to derive 
the empirical Bath's law 111 which states that the average difference in magnitude between 
a mainshock and its largest aftershock is 1.2 regardless of the mainshock magnitude, within a 
completely consistent theory taking into account all generations of triggered events. Our present 
results significantly improve on the numerical results of 12] by demonstrating the essential roles 
played by the cascade of triggered events and the proximity to criticality in order to obtain Bath's 
law and by providing new improved constraints on the key parameters of the ETAS model. In 
addition, we extend Bath's law which is a statement on the average magnitude difference between 
the mainshock and its largest aftershock by giving the full distribution. Our theoretical framework 
also allows us to calculate precisely the probability that the largest aftershock turns out to be larger 
than its source, a situation which is usually interpreted as the source and all events before the largest 
aftershock being its foreshocks, the largest aftershock being re-interpreted as the mainshock of the 
seismic series. 

The paper is organized as follow. The next section 2 recalls the definition of the branching 
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model of triggered seismicity. Section 3 presents the generating probability function (GPF) and 
results on the statistics of the largest aftershock among aftershocks of the first generation. The 
GPF for first-generation aftershocks is generalized to aftershocks of all generations in section 4. 
This allows us to predict that the distribution of magnitudes of the largest aftershock over all 
aftershock generations is renormalized in the critical regime. This renormalization provides a novel 
way to calibrate the productivity parameter. Section 5 puts together previous results to calculate 
the average difference in magnitude between the mainshock and its largest aftershock over all 
generation. In the critical regime, Bath's law is shown to hold. The value of the average difference 
in magnitude allows us to offer new improved constraints on the two key parameters of the ETAS 
model, the critical branching ratio n and the productivity exponent a. Section 6 concludes. 



II. THE EPIDEMIC-TYPE AFTERSHOCK SEQUENCE (ETAS) BRANCHING 
MODEL OF EARTHQUAKES 

Consider an earthquake of magnitude p, which we refer to as a mainshock to mean that we 
are interested in the earthquakes that it triggers (aftershocks). According to the ETAS model, it 
generates a random number i?^ of first generation aftershocks, which has Poissonian statistics, 

characterized by the conditional average number 

N, = kM. Mp) = 10"("-"^°^ (2) 

Here itlq is the minimum magnitude of earthquake capable of triggering other earthquakes, and k 
is a constant. The expression Q for yu(p) is chosen in such a way that it reproduces the empirical 
dependence of the average number of aftershocks triggered directly by an earthquake of magnitude 
m (see jisl \i\ and references therein). Expression Q gives the so-called productivity law of a 
given mother function of its magnitude p. 

The ETAS model requires the specification of the Gutenberg-Richter (GR) density distribution 
of earthquake magnitudes 

p{m) = b In(lO) lO-^('"-™o) , m ^ mo , (3) 
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such that p{x)dx gives the probabihty that an earthquake has a magnitude equal to or larger 
than m. This magnitude distribution p{m) is assumed to be independent of the magnitude of the 
triggering earthquake, i.e., a large earthquake can be triggered by a smaller one j^. The cumulative 
{V{m)) and complementary cumulative (Q(m)) distributions corresponding to (jS)) are 

V{m) = 1 - Q(m), Q(m) = io-^("o) = [^^{m)]-^ , 7 = 6/« • (4) 

The ETAS model is defined by the conditional Poisson intensity given the average rate of 
seismicity at time t and position r conditioned on all past earthquakes 

A(t, r) = s{t, r)+ ^(^ - - ^i) , (5) 

i I ti<.t 

where s{t,r) is the average Poisson rate of spontaneous earthquake sources ("immigrants" in the 
language of epidemic branching processes) at position r and at time t. The sum is over all past 
earthquakes: each earthquake is characterized by its occurrence time ti, its magnitude and its 
location rj in the catalog. The two kernels — ti) and 0(r — rj), whose integrals with respect 
to time and space respectively are normalized to 1, describe the contribution of the earthquake at 
time ti and position to the seismic intensity at time t in the future and at position r. 

III. STATISTICS OF THE LARGEST AFTERSHOCK AMONG AFTERSHOCKS OF 
THE FIRST GENERATION 

A. Generating probability function (GPF) of aftershocks of the first generation triggered 
by an arbitrary mainshock 

The Poissonian statistics of the aftershocks of first generation implies that the generating prob- 
ability function (GPF) of their numbers reads 

00 

e^{z\p) = {z'^p) = J2 Ppij-)z'' = e"^(^)("-^) , (6) 

r=0 

where we have used (0) and the angle brackets correspond to taking the statistical averaging. 

Let Ml, M2, . . . Mfii be the random magnitudes of the aftershocks of first generation triggered 
by the mainshock. Let us consider the statistical average defined by 

ei{z,m\p) = {z''pl[H{m-Mk)) , (7) 

k=l 
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where H is the Heaviside function. Using the Poissonian statistics of the random number of 
aftershock numbers and the Gutenberg-Richter law for their magnitudes, we obtain 

oo 

ei(z, m\p) = ^ Ppir)z'[V{m)Y = g-'^t^^M-il (8) 

r=0 

Notice that Qi{z,m\p) can be rewritten as 

ei{z, m\p) = {z^p^'^^Him - M])) , (9) 

where is the largest magnitude over all aftershocks of the first generation triggered by the 
mainshock and Rp{m) is the number of aftershocks for those realizations of aftershocks in which 
no aftershocks of the first generation exceed the magnitude m. 

The interest in Qi{z,m\p) in (j8l9p lies in particular in the fact that, for z = 1, it reduces to 
the probability Pi{m\p) that the largest magnitude among all aftershocks of the first generation is 
smaller than m: 

Pi{m\p) = Giiz = l,m\p) = Pr {M^ < m} = e-'^'^('')2(m) ^ ^^g) 

B. GPF of aftershocks of the first generation triggered by a spontaneous source 

Consider now some spontaneous source (contributing to the term s{t, r) in (0)), with magnitude 
Mq. According to the ETAS model, it triggers its own aftershocks sequence independently of all 
other sequences. Let 

Mo, Ml, M2,...M^i (11) 

be the random sequence of magnitudes (including Mq) of the first generation aftershocks triggered 
by the spontaneous source. Then, analogously to (jH)), the GPF 

Qi{z, m) = {z^'^"''>H{m - M^)) , (12) 

describing the statistics of the number R^{m) of aftershocks of first generation (triggered by the 
spontaneous source) and their largest magnitude M^ among the list (jllj) . is equal to 

ei{z,m) = F[l - zV{m)] - [/i(m)]"^F(^(m)[l - zV{m)]) , (13) 

where 

F{y)=^K^y^T{-^,Ky) . (14) 
5 



For m oo, the GPF Qi{z, m) given by (fT^ reduces to the standard GPF of the random number 
Ri of aftershocks of the first generation triggered by some spontaneous source, which reads 

e^{z) = ei{z,m = oo) = F{l-z) . (15) 

In the following analysis, the function F{y) in ()14|) plays a crucial role. It is thus useful to 
state some of its analytical properties. For 1 < 7 < 1.5 and y < 0.2, it can be represented rather 
accurately by 

F{y):^l-ny + (3y\ /3 = - (^n^^ r(l - 7) . (16) 

In the theory of aftershocks branching processes, the branching ratio, equal to the average 
number of first generation aftershocks triggered by some mother earthquake, plays a fundamental 
role since it controls the subcriticality versus supercriticality of the process. It is defined as 



dQiiz) 

n = 



-, — a 
2=1 



(17) 



dz 

Thus, for given 6, a and n, the constant k in relations (jH)) and (jl(J|) can be replaced by 

K = K(7,n)=n(^l-|) . (18) 

Taking z = 1 in (|13|) obtains the cumulative distribution function (CDF) Pi{m) of the largest 
magnitude among the sequence (fTT|) of aftershocks including their spontaneous source 

Pi = F{Q)-fi-^F{fiQ). (19) 

The GPF 9i(z,m) in (fT^ describes the statistics of the spontaneous source and its first gen- 
eration aftershocks, such that all magnitudes, including the magnitude of the spontaneous source, 
are smaller than m. In the following, we will consider the possibility that the largest aftershock 
may be larger than the source, a situation which is known in the seismological literature as the 
occurrence of foreshocks (see j^Q,^^ and references therein). The corresponding GPF, averaged 
over all possible foreshock's magnitudes, is denoted as 6i{z,m) and is obtained formally from (jl3|) 
by taking the limit /i — > 00: 

9i{z,m) = F[l-zV{m)] . (20) 
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C. Magnitude of the largest aftershock of the first generation aftershocks 



Before analyzing the conditions under which the empirical Bath's law can be obtained from the 
ETAS model, it is useful to ask what is its analog when restricting the set of aftershocks to the 
first generation of events triggered by the source. This will provide a reference point against which 
to gauge the impact of the multiple generations of aftershocks on Bath's law. 

We start from expression (ITUI) giving the CDF Pi(m|p) of the largest magnitude among all 
aftershocks of the first generation. Substituting ()2l4p in (fTIH). we obtain 

Pi{m\p) = G[wo{m - p + vo)] , (21) 

where 

vo{p) = (l - (p - mo) + ^ logio ( ^(^ \ ) > ^i;o = &lnlO, (22) 

and 

G{x) = exp (-e"^') (23) 

is the well-known limiting extremal Gumbel CDF. 

It follows from (PT|) in particular that the probability density function (PDF) of the difference 

a; = p - m; (24) 

between the source (mainshock) magnitude and the magnitude of its largest aftershock of the first 
generation is equal to 

mp) = Wo g[woivo - 6)] , (25) 

where 

g{x) = exp(— X — e~^) (26) 

is the PDF associated with the CDF ()23|) . Note that the shape and variance of the PDF 
does not depend on the mainshock magnitude p. Only its mode fo(p) (most probable value of the 
difference ()24j) ) depends on p and increases linearly with it according to (j22I)- 

These results treat all aftershock sequences on the same footing and in particular include se- 
quences which have zero aftershocks. In a real data analysis, the statistical properties of the 
difference ()24p are obtained conditioned on the observation of at least one aftershock, which re- 
quires a modification of the expressions above. We are interested in modifying the CDF (fTO|) to 



eliminate the cases where Rp = 0. This corresponds to obtaining the CDF of the largest magni- 
tude of first generation aftershocks under the condition that the mainshock triggers at least one 
aftershock: 

I _ g-KM(p)S(m) 

Pi(m|p;l) = l-gi(m|p;l), Q,(m\p;l) = . (27) 

The corresponding PDF of the difference (j24p reads 

/•W=»»T=T^- — <^</'-"o. (28) 
The conditional CDF ()27|1 differs significantly from the unconditional one only if the probability 
g-K/i(p) ^j-^at there are no aftershocks is close to 1. This occurs for small mainshock magnitudes. 
How small should the mainshocks be for this difference to be important? Let us define a magnitude 
threshold po by 

Kp(p)~2 ^ pQ = mo + - logio ( N ) • (29) 

a \n{0 — a) J 

For mainshock magnitudes p > po, the conditional CDF ()27p does not differ significantly from 
the unconditional one (jlUj) . In this case, relations (j25|) and (|28p are approximately equal, and the 
looked for distribution of the difference ()24j) can be taken to be the PDF ()25j) where 6 G {—oo, oo). 
Thus, for p > Po, the average of the magnitude difference can be approximated by 

/oo 
6 g[wo{vo - 6)]d6 = Vo{p) - u/wo , (30) 
'OO 

where p ~ 0.5772 is the Euler constant. Figure 1 shows the exact average difference A^m calculated 
with (j28p . its approximation (j30|) valid for sufficiently large mainshocks p > po and the most 
probable value A^m,, = vo{p) of the difference in magnitude between the mainshock magnitude p 
and its largest aftershock. This figure is typical of the strong dependence found for all reasonable 
values of the parameters and distinguishes this result from the empirical Bath's law (which gives 
a constant value independent of p). 



IV. STATISTICS OF THE LARGEST AFTERSHOCK AMONG AFTERSHOCKS OF 
ALL GENERATIONS 

A. GPF of the aftershocks over all generations triggered by a spontaneous source 

Due to the mutual statistical independence of different branches of triggered earthquakes in the 
ETAS model, one can easily generalize the results for the largest aftershock of the first generation 
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to derive the statistical properties of the largest aftershock over all generations. 

Within the ETAS branching model, taking into account all aftershocks of all generations which 
are triggered by the mainshock amounts to replacing in the r.h.s. of equation (jH)) the PDF V{m) 
of the magnitudes of single aftershocks by the GPF 0(z, m) for all aftershocks triggered by some 
spontaneous source which have (together with the source) magnitudes smaller than m. As a result, 
we obtain the sought GPF of the number of all aftershocks triggered by the mainshock conditioned 
on all magnitudes to be less than m as 

Q{z, m\p) = H{m - Mp)) = e''Mp)[-e{.,m)-i] _ ^3^^, 

Here, Mp is the magnitude of the largest aftershock. In particular, the complementary CDF of the 
magnitude of the largest aftershock reads 

Q{m\p) = {H{Mp -m)) = l- e-'^^(^)Q('^) . (32) 

Similarly, the functional equation for the GPF Q{z, m) is obtained by replacing in p3|l both 0i 
and V{m) by Q{z,m): 

e{z, m) = F{1 - ze{z, m)) - fi-^{m)F[ii{m){l - zB^z, m))] . (33) 

For -2 = 1, equation (j33|) reduces to an equation for the CDF P(m) of the magnitude M of the 
largest event (including all aftershocks and the source) defined as 

e{z = l,m) = P{m) = Ft{M <m} . (34) 

This equation reads 

P = FiQ)-fi-^F{fiQ), Q{m) = 1- P{m) . (35) 
Notice that in the limit m — > 00, Eq. (|33|) reduces to the well-known functional equation 

9 = 91(^9) (36) 

for the standard GPF 

e(z, m = 00) = &{z) = (z^) (37) 
of the random number R of all aftershocks triggered by some ancestor, which has been studied in 



Similarly to the reasoning leading to (P7j) . the conditional probability that the magnitude of 
the largest aftershock exceeds m, under the condition that the mainshock triggers at least one 
aftershock, is 

1 _ g-K/i(p)<9(»n-) 

Q[m\p-l) = —^-—^^, Q{m) = l-P{m). (38) 

It is also of interest to obtain the GPF 0(-2, m) of the number of aftershocks of all generations 
with magnitudes smaller than m which are triggered by some spontaneous source of arbitrary 
magnitude. It is given by replacing in the r.h.s. of Eq. (j2(J|) V{m) by the GPF Q{z,m) given by 
(jHHj) which yields 

e{z,m) = F[l- ze{z,m)] . (39) 

For z = 1, this gives the probability P{m) that the magnitude of the largest aftershock triggered 
by an arbitrary spontaneous source is smaller than m: 

P(m) = F[Q{m)] . (40) 

B. Distribution of the magnitude of the largest aftershock of a spontaneous earthquake 
source of arbitrary magnitude 

All quantities defined above require the knowledge of Q(m), which has a straightforward sta- 
tistical meaning: it is the complementary CDF of the magnitude of the largest aftershock of a 
spontaneous source of arbitrary magnitude (in other words, over all possible source magnitudes). 
It is easy to calculate Q{m) by solving equation numerically. We can also use the algebraic 
approximation (fT?)|l of the function F{y) to obtain an explicit and rather precise analytic expression 
of Q{m). Indeed, substituting (fTBj) into (P^j) obtains 

- T. / \ ■ (41) 

^ ' (1 - n)[/i(m)]T + n/i(m) ^ ^ 

This expression shows that there is a cross-over magnitude mc, given by 

/ X / ^ \ 1/(7-1) 1 , / n \ , ^ 

/i mc ~ ^ mc ~ mo + logio 1 , 42 

separating two regimes with different power laws for Q{m). The first regime 

Q(m) ~ lio-°(™-™o)^ m < mc , (43) 
n 
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corresponds to a complementary CDF decaying slower than the Gutenberg-Richter law ()3l4j) . for 
a < b. In the critical case n = 1, = oo and this regime PHj) holds for any m > ttiq. The second 
regime recovers the Gutenberg-Richter law 

g(m) ~ _J_iQ-b{m-mo) ^ m>m^. (44) 
1 — n 

Figure 2 shows the logarithm (in base 10) of the complementary CDF Q{m) as a function of m — 
and the two power law asymptotics ()42|) and PHj) . We have thus shown that the Gutenberg-Richter 
law can be renormalized from a bare exponent 6 to a smaller exponent a when the distribution 
is restricted to the set of largest aftershocks of spontaneous earthquakes of arbitrary magnitudes. 
This renormalization of the 6-value from 6 to a is intrinsically a cascade phenomenon. In other 
words, it results from the existence of a cascade of triggered earthquakes over many generations 
as shown in details in ^. This renormalization proceeds W a mechanism similar to that of the 
p-value of the Omori law from a value 1 + to 1 — ^ It is different from the mechanism 

leading to an exponent b — a for the asymptotic branch of the Gutenberg-Richter distribution of 
all foreshocks [la]. 

This prediction ()4Hj) offers a novel method for measuring the key exponent a controlling the 
productivity or triggering efficiency of earthquakes as a function of their magnitude, according 
to (j21). What is needed to implement this new method is a declustering technique to identify the 
spontaneous sources and their largest aftershocks. The statistical declustering technique of Zhuang 
et al. [3| seems to be particularly suitable for this purpose. 

Due to the independence between the aftershock sequences of different sources in the ETAS 
model, it is straightforward to obtain the probability distribution of the largest triggered events 
among a set of r spontaneous sources. The corresponding complementary CDF, giving the prob- 
ability that the largest event triggered over all generations by r sources is larger than m, is equal 
to 

Q{m\r) = 1- F''[Q{m)] . (45) 

Figure 3 plots Q{m\r) as a function of m — mo, for r running from 1 to 15 for a = 0.8, 6=1 
and for two values of n, n = 0.9 and n = 1. It is clear that most of the largest triggered events 
in aftershock sequences are very small, simply due to the interplay of two factors: most random 
sources are themselves small and have a small productivity, and the Gutenberg-Richter distribution 
makes it much more probable that all triggered events have small magnitudes. 
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C. Distribution of the magnitude of the largest aftershock of a spontaneous earthquake 
source of fixed magnitude p 



Rather than considering arbitrary source magnitudes, it is interesting to determine the comple- 
mentary CDF Q{m\p,r) of the magnitude of the largest event triggered by r spontaneous sources 
with fixed magnitudes pi, p2, . . . , Pr- It is closely approximated by 



Q{m\p, r) ~ 1 — exp 



i=l 



(46) 



In the following, we restrict our analysis to the case of a single r = 1 spontaneous source which 
fixed magnitude p. In this case, expression (|46|) transforms into (j32j) . 

For m < rric, where the cross-over magnitude rric is defined by (j42j) . expression (|32j) can be 
simplified by replacing the exact Q{m) by the approximation which gives 

Q{m\p) ~ 1 — G {wi{m — p + Vi)) rriQ < m < rric . (47) 

where 

vi = -\og^A- ), M;i = alnlO. (48) 

Expression ()47|1 can be further simplified into 

Q(m\p) ^ - 10-"("-^) = (k 10°(^-™o)) X (- io-"("o)) , (49) 
n ^ '' \n J 

in the tail of Q{m\p\ i.e., for p — Mog^g f ^ '^^^ re-writing of Q{m\p^ under the form shown 
by the last equality in clarifies its origin: the first factor k 10"^''"™'°) is nothing but the 
productivity law the second factor ^ ]^Q-a(m-nio) jg ^j-^g renormalized Gutenberg- Richter law 

For m > VTif.^ we obtain another approximation for Q(m|p) by replacing in the r.h.s. of equation 
the complementary CDF Q{m) by the approximation ()44|) . which yields 



Q{m\p) 1 - G {w2[m - p + V2{p)]) m>mc, (50) 

where 

..= (l-^)(p-™„) + ilog,„(;|-^) . ». = MulO. (51) 
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The following approximation including both regimes m < rric, m > rric and laws 1)471501) is 
obtained from (jH^ by replacing Q{m) by the approximation ()41|) : 



Figures 4 and 5 present the dependence of the complementary CDF Q{m\p) as a function of the 
magnitude of the largest event triggered by a spontaneous source of fixed magnitude p, for four 
different values of p. The figures show the exact Q(m|p) obtained numerically, its approximation 
(jH2I) (which is actually undistinguishable from the exact one) and the universal approximation 
(^Tj). The comparison between Figure 4 (for a = 0.8) and Figure 5 (for a = 0.9) shows that the 
approximation (|47j) becomes more and more precise as a becomes closer to b. 

As a bonus, we obtain the probability Q*{p) that the magnitude m of the largest aftershock 
exceeds the source magnitude p, a situation which is usually classified in seismic catalogs by saying 
that the largest triggered event is the mainshock and the spontaneous source that initiated the 
sequence and all triggered events before the largest aftershock are foreshocks. Indeed, Q*{p) is 
nothing but 



Interestingly, in the regime mo < p < ttLc (i.e., for n sufficiently close to 1 and/or a close to b) for 
which (jTTj) holds, we obtain 



which is independent of p and of the branching ratio n. This approximation is all the better, the 
closer a is to b. Figure 6 shows the exact Q^.{p) as a function of p, which can be compared with 
the constant (jHH). One can observe that, at least for a = 0.95, Q*{p) is actually quite close to 
the constant (jK^ over all possible magnitudes p > niQ. We propose to call the prediction (jK^ 
the "second Bath's law for foreshocks: the probability that a main earthquake turns out to be the 
foreshock does not depend on its magnitude p (more generally, the distribution of the difference 
p — m does not depend on p). 




(52) 



Q.{p) = Q{p\p-A) ■ 



(53) 




(54) 
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V. DERIVATION OF BATH'S LAW 



The derivation of the distribution (jHHj) and the approximations (jTTj) and ()50|) allow us to derive 
Bath's law by calculating the statistical average 

Apm = p - {Mp) , (55) 

where Mp denotes the magnitude of the largest aftershock among all events of all generations 
triggered by the source of fixed magnitude p. Recall that Bath's law states that Apm is independent 
of p and equal to 1.2. 

Within the ETAS model, the exact value of ApUi is obtained as 



oo 



ApTrii ~ f 1 = B = — 

Wi a 



log 



10 



(57) 



ApTn = p — Q{m\p; l)dm (56) 

J mo 

where Q{m\p;l) is given by (jHH|) and the r.h.s. of (jKHjl expresses the fact that the average is 

performed over sequences with at least one aftershock. 

The two regimes m < rric giving the asymptotic (P7j) and m > rric giving the asymptotic ()5Up 

provides two asymptotic expressions for Apm. Indeed, calculating (Mp) using the approximation 

(j47j) and neglecting the boundary effects (i.e., supposing that m G {—oo, oo)) obtains 

' b \ _ V ' 
b-a) In 10 

which is independent of p. Recall that u ~ 0.5772 is Euler constant. In the following, we call B 
defined in fj57|) the Bath's constant. Note that this regime m < rric corresponds to the critical 
branching regime of n close to 1 (for a fixed magnitude p) and expresses the full effect of the 
cascade of triggered events over all possible generations. It is remarkable that the theory of the 
ETAS branching model predicts the first part of Bath's law that the average of the difference 
between the magnitude of a mainshock and its largest aftershock is independent of the mainshock 
magnitude. The specific value of Bath's constant B depends only two parameters, the 6-value of 
Gutenberg- Richter distribution and the productivity exponent a. 

The second asymptotic for m > rric corresponds to using (p])) to estimate (Mp), which yields 



Apma ~ t;2 - — = f 1 - ^) (p - mo) + 



log 



^1 — n)b\ V 
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n{b-a)l In 10 



(5^ 



W2 \ b J b 
Note that ApUi^ is increasing with p as in expression (jHUj) corresponding to taking only into account 
aftershocks of the first generation. This is natural since the asymptotic for m > rric corresponds to 
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n relatively far from 1 (for a fixed p), i.e., far from the critical branching regime, such that only a 
few generations play a significant role in the population of aftershocks. This asymptotic ()58|) is also 
identical to the expression (5) of derived by using the statistical average of the total number 
Ai'aft of aftershocks of all generations triggered by a source of fixed given magnitude. Thus, the 
difference between this approximation ()58|) (and expression (5) of 1^) and the exact expression 
fl56|) and its critical universal asymptotic (jFfj) can be traced back to the difference between the 
following two kinds of averages: (ln[A''aft]) and In(A^aft)- 

Figures 7 and 8 show the exact average magnitude difference (jSUj) as a function of the mainshock 
magnitude p for 6 = 1 and different values of the branching ratio n, for a = 0.9 (Figure 7) and 
a = 0.95 (Figure 8). As expected from the condition m < rric or n closer to 1 at fixed p so that 
rric is all the larger according to (021), the largest values of n give almost constant values of Apm 
in agreement with the prediction (jFf|) . For smaller n's, we can observe a slow cross-over to the 
second asymptotic (j3H|) . By comparison between Figure 7 and Figure 8, it is clear that low values 
of a are not compatible with Bath's law. As confirmed with similar figures obtained for smaller 
a's, a value of a at least equal to 0.9 seems necessary to obtain a dependence of Apin roughly 
independent of p over a large magnitude range. This bound is compatible with some previous 
studies |l4, Q but in disagreement with others 18, 20 1. But, we should remark that such 
heterogeneity in reported values of key parameters of the ETAS model such as the productivity 
exponent a could be due to the bias resulting from imperfect account of unobserved seismicity 
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. Figure 



below the completeness threshold, which may play a dominant role as explained in 
8 also shows that, the closer is to 1, the more independent is Apm with respect to the mainshock 
magnitude p. But, due to inherent fluctuations in empirical data, n can be as low as n = 0.8 for 
a ~ 0.95 and ApUi would still be slowly growing between 1.1 and 1.3 over a large magnitude range 
of p, so that Bath's law would be approximately verified. 



VI. CONCLUDING REMARKS 



Using the ETAS branching model of triggered seismicity, we have shown how to calculate exactly 
the average difference between the magnitude of a mainshock and the magnitude of its largest 
aftershock over all generations. This average magnitude difference is found empirically to be 
independent of the mainshock magnitude and equal to 1.2, a universal behavior known as Bath's 
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law. We have developed the mathematical formulation in terms of generating probability functions 
which allow us to obtain exact equations and useful approximations to understand the physical 
basis for Bath's law. In particular, we find that the constancy of the average magnitude difference 
(to a value that we term Bath's constant) is associated with the critical regime of the ETAS 
branching process. Allowing for error bars ±0.1 for Bath's constant value around 1.2, our exact 
analytical treatment of Bath's law provides a new constraint on two key parameters of the ETAS 
model, namely the productivity exponent a and the branching ratio n: a > 0.9 and n > 0.8. 
We have suggested a novel method for measuring a based on the predicted renormalization of the 
Gutenberg-Richter distribution of the magnitudes of the largest aftershock. To implement this 
method, statistical declustering techniques can be used to identify the spontaneous sources and 
their largest aftershocks. We have also proposed the "second Bath's law for foreshocks" that the 
probability that a main earthquake turns out to be the foreshock does not depend on its magnitude 
P- 
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Fig. 1: Exact average of the difference Aim = p — (Mp) obtained using H28() (bottom curve 
bending down for small p — tuq), its large magnitude approximation (|3U|) (bottom straight 
line) and the difference A^m* = vq{p) between the mainshock magnitude p and the mode of 
the magnitude of the largest aftershock among all aftershocks of the first generation (upper 
straight line), for n = 0.9, a = 0.8 and 6 = 1. 
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Fig. 2: Plot of the decimal logarithm of the exact CDF Q{m) and its approximation (|41l) 
(which actually coincide). Straight lines correspond to the asymptotics and 
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0.5 1 1.5 2 2.5 3 3.5 4 

Fig. 3: Complementary CDF (|45j) of the magnitude of the largest event triggered by r 
spontaneous sources with random magnitudes chosen according to the Gutenberg-Richter 
distribution. The different curves correspond to r = 1; 3; 5; 7; 9; 11; 13; 15 from bottom to top. 
Each couple of curves corresponds to n = 0.9 and n = 1 respectively. 
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Fig. 4: Exact complementary CDF's Q{m\p) given by 1)321) (lower curves), their approxima- 
tions ()52() (which are actually undistinguishable from the exact functions) and the universal 
approximations (|47|) (upper curves), for n = 0.9, 6 = 1 and a = 0.8 (7 = b/a = 1.25) for four 
different values of the spontaneous source magnitude p. 
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Fig. 5: Same as Figure 4 for a = 0.9 (7 = 1.11). 
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Fig. 7: Plot of ApTn given by the exact expression ()56l) for a = 0.9, 6=1 and for different 
branching ratio. Up to down n = 0.8; 0.85; 0.9; 0.95; 0.99. With n tending n to 1, the average 
magnitude difference become closer to the theoretical Bath's constant B = 0.83. 
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Fig. 8: Same as Figure 7 (except for the magnification) for a = 0.95 (7 ~ 1.05) giving 
B = 1.11. 
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